Two-field Models of Dark Energy with Equation of State Across -1 
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In this paper, we study the possibihty of building two-field models of dark energy with equation 
of state across -1. Specifically we will consider two classes of models: one consists of two scalar 
fields (Quintessence + Phantom) and another includes one scalar (Phantom) and one spinor field 
(Neutrino). Our studies indicate to some extent that two-field models give rise to a simple realization 
of the dynamical dark energy model with the equation of state across w = —1. 
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I. INTRODUCTION 

Recent observational data Q, |^ strongly indicate 
that the Universe is spatially flat and accelerating at 
the present time. In the frame of Friedmann-Robertson- 
Walker (FRW) cosmology, the acceleration may be at- 
tributed to some mysterious source called dark energy. 
The simplest candidate for dark energy seems to be a 
small positive cosmological constant, but it suffers from 
difficulties associated with the fine tuning and coinci- 
dence problem. An alternative is a dynamical scalar field, 
such as Quintessence ^] or Phantom The Phantom 
field violates the energy conditions, which leads to many 
interesting cosmological phenomena I3| . 

Despite the current theoretical ambiguity for the na- 
ture of dark energy, the prosperous observational data 
(e.g. supernova, CMB and large scale structure data and 
so on ) have opened a robust window for testing the re- 
cent and even early behavior of dark energy using some 
parameterizations for its equation of state. The recent 
fits to current supernova(SN) la data, CMB and LSSQ 
find that even though the behavior of dark energy is to 
great extent in consistency with a cosmological constant, 
an evolving dark energy with the equation of state w 
larger than -1 in the recent past but less than -1 today 
is weakly favored. If such a result holds on with the ac- 
cumulation of observational data, this would be a great 
challenge to the current cosmology. 

The evolving dark energy with an equation of state w 
crossing -1 during its evolution was firstly advocated and 
named as Quintom in Ref. The Quintom models of 
dark energy arc different from the Quintessence or Phan- 
tom in the determination of the evolution and fate of the 
universe. Generically speaking, the Phantom model has 
to be more fine tuned in the early epochs to serve as dark 
energy today, since its energy density increases with ex- 
pansion of the universe. Meanwhile the Quintom model 
can also preserve the tracking behavior of Quintessence, 
where less fine tuning is needed. The Quintom model 
with an oscillating equation of state considered in Qcan 
lead to the oscillations of the Hubble constant and a re- 
curring universe, which in some sense unifies the early 
(Phantom) inflation ^3 and current acceleration of the 
universe. This oscillating Quintom would not lead to a 
big crunch nor big rip. The scale factor keeps increas- 



ing from one period to another and leads naturally to 
a highly flat universe. Since in this model the universe 
recurs itself and we are only staying among one of the 
epochs, the coincidence problem in some sense is recon- 
ciled. Ref.TH has considered a variation of this oscillat- 
ing Quintom model and studied its constraints from SN, 
CMB and LSS. 

There have been some efforts made on the model build- 
ing of the Quintom dark energy based on field theory. 
First of all, a single scalar field with the canonical kinetic 
term (Quintessence) or a negative kinetic term (Phan- 
tom) will not be able to give rise to the equation of state 
across -1. The similar conclusion has also been obtained 
for the k-essence models ^31. Beyond the single scalar 
field theory, Refs. 8] and [l^ have proposed an explicit 
model of Quintom with two scalar fields, one being the 
Quintessence and another being the Phantom field. This 
type of model can easily lead to a scenario where at early 
time the Quintessence dominates with w > —1 and lately 
the Phantom dominates with w < —1, satisfying current 
observations. Some recent relevant studies are given in 

The Quintom model considered in refs. |1] and |0| does 
not include the interaction between the two scalars and 
is not the most general one. In this paper, we revisit 
this type of two-field Quintom model by introducing a 
interaction term. We will study in detail the cosmo- 
logical evolution of this model. In ref.|^3 it has been 
shown that in the two-field Quintom model the scaling 
solution dominated by the Phantom field will be a late- 
time attractor in the absence of the interactions. In this 
paper by a explicit calculation we will show the inter- 
actions do not affect the Phantom-domination attractor 
behavior. In addition our studies also show that this 
class of models will provide some interesting possibilities 
of the evolution of the equation of state which have not 
been considered in refs.^ and jl3j|. Furthermore we will 
propose in this paper a two-field Quintom model with 
a scalar and a fermion which specifically we take to be 
the neutrinos and study its evolution of the equation of 
state. This paper is organized as follows: in section II we 
study in detail the model with one Quintessence and one 
Phantom; in section HI we present our model of Quin- 
tom with a Phantom and neutrinos; section IV is our 
conclusion and discussions. 
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II. QUINTOM MODEL WITH PHANTOM 
FIELD AND QUINTESSENCE FIELD 

In this section, we study a Quintom model with two 
scalar fields, in which one is the Phantom field and an- 
other is the Quintessence field. The Lagrangian for such 
a coupled system is given by 



^ — 1^4, + ^cr + Vint , 



where 



C4 = -T:d^4>d^(t>~V{4>), 



(1) 



(2) 



with V{(j)) being the potential for the Phantom field and 



(3) 



with V{(j) the potential for the Quintessence field, and 
Vint denotes the interaction between the Phantom and 
the Quintessence fields. The equation of state of this 
system is 



w 



-\(f' + \a''-V{<j))-V{a)-V^nt 
-\4>'' + h^^ + V{^) + V{a) + V,nt 



(4) 



The evolution equations of the fields and the fiuid for 
a spatially fiat FRW model are 



cr + ZHa H ; ! — = 0, 

da da 



(5) 

(6) 
(7) 



where is the density of fluid with a barotropic equation 
of state = (7 — where 7 is a constant in the range 
of < 7 < 2. 

Now we study the cosmological evolution of the model 
with lagrangian (1). The authors of ref. [3 have con- 

_Xi-i— 

sidered a model where V(4>) = V^qc ^^p' , V{a) = 
-A2 " 

Vao e ^'p' and show that the late time behavior of this 
model is a Phantom-domination attractor. Here we in- 
troduce an interaction term Vint — ^{V{4>)V{a)Y^''^ and 
study its effect on the cosmological evolution. We follow 
closely the conventional phase-plane analysis for the spa- 
tially flat FRW models in [isl and the detailed studies 
on models with multi coupled Quintessence fields in Ref. 
[T^ . Defining firstly the variables 



a 



(8) 



the evolution Eqs. (5-7) can be written to an autonomous 
system 

%/6, 



-3x^(1 + xl-xl-^)- ^Xivl 



V6 



XiXy<pya, 



(9) 



^y<phxl + xl + '^-) - ^Xix^y^, (10) 



V6, 



o n _L 2 2 , \/6 

—.iXa(i-\-x^ — x^ 

+—X2Xy4,ya, 

2 

o I 2 , 2 I 7^ \ 
iy<y(.'X^ + x^ + —) 

or 2,2, 7^^ \ 3 



(11) 



V6 



X2XaVa: (12) 

(13) 



where a prime denotes a derivative with respect to the 
logarithm of the scale factor. The critical points corre- 
spond to the fixed points where x'^ = 0, x'„ = 0, y'^ = 0, 
y'^ = 0, z' — 0, which have been calculated and given in 
Table 1. 
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FIG. 1: Plot of the equation of state w as a function of the 
scale factor Ina. 



Eqs. (9-13) can be reduced to four independent equa- 
tions. To study the stability of the critical points, we sub- 
stitute the linear perturbations about the critical points 
into these independent equations and keep terms to the 
first-order in the perturbations. The four perturbation 
equations give rise to four eigenvalues. The stability re- 
quires the real part of all eigenvalues be negative (see Ta- 
ble II for the eigenvalues of perturbation equations and 
the stability of critical points). 

From Table I and II, one can see that even if there ex- 
ists the interaction between the two fields, the Phantom- 
dominated solution is still a late-time stable attrac- 
tor. If this coupled system is initially dominated by 
the Quintessence field, it will eventually evolve into the 
Phantom-domination phase. Thus an equation of state 
across —1 will be inevitable. This provides a natural 
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Label 
T. 

2. 

3. 
4. 

5. 



X4, 
-7 



y<t> 



Xa, — 1 



A 








A2_ 



37 








37(2-7) 



1 37 



TABLE L The list of the critical points. 



Label 



mi 



m2 



rn-i 



mi 



3(1- ^Aii^,) 
-i(6 + Af) 

-3(1- f) 

37 
2 

-(3-^) 



3(1 - ^Mx^^) 
-37 - A? 
-37 + Ai 
-3(1 -i) 



Stability 



1. 

2. 

3. 
4. 

5. 



-6(1 - ^)x1^ 
_^ 

-'2 

ii 
2 

2 



6(l-i)<c 
-1(6 + A|) 

-3(1- f) 
-3(1 -i) 



-3(2-7) 



87A^ 

Ta^ 



-247^ 



unstable 

stable 

unstable 
unstable 



-3(2-7) 



87Aj- 



) unstable 



TABLE IL The eigenvalues and stability of the critical points. 



realization of the Quintom scenario. As an illustra- 
tion we in Fig. 1 plot the evolution of the equation 
of state w as a function of the scale factor where in 
the numerical calculation we have taken V^^ = V^^, = 
0.35 X 10-46(Gel^)4, Ai = A2 = A = 1. 

Before concluding this section we point out some in- 
teresting behavior of the equation of the state of the two- 
field models. In Fig. 2 and Fig. 3, by choosing potentials 
and specific values of the model parameters we show that 
the equation of state could be oscillated. And the oscil- 
lations mainly occur in Fig. 2 in the late time when 
redshift 2; < 0, while Fig. 3 give rise to one example more 
interesting, in which the oscillations across —1 occur in 
the near past, which might lead to some effects testable 
observationally. 
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FIG. 2: The state equation «; as a function of scale factor In a 
for the potential V = \\ cos (Ci^) + -^2 cos(Q-j^|-^) + \(t>'^a^ , 

where ^ = 0.5, a = 1, Ai = A2 = 2.47 x IQ'^^GeV*. The four 
lines from the top to the bottom at Ina = correspond to 
A = (0.2,0.3,0.4,0.5) x 10"^^° respectively. 
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FIG. 3: The state equation w as a function of scale factor 
Ina for the potential V = ^mlcf)^ + ^mla'^ + XcfPa'^, where 
\ml = 5.9 X 10"*^GeV^ \mi = 1.19 x lO-^^GeV^ and A = 
50 X 10"^^°. 



III. QUINTOM MODEL WITH PHANTOM 
FIELD AND NEUTRINO 

The two-field models of Quintom dark energy consid- 
ered above consist of two scalar fields. In the follow- 
ing we study a model where the Quintessence in |^ is 
replaced by the neutrinos. In this model the neutrino 
will become a part of the dark energy. There have been 
a lot of studies recently in the literature on the cou- 
pled system of neutrinos interacting with the dark energy 
scalars|i3 EHIlllSIIllllllllMllI, however in all 
of these studies the scalar fields have canonical kinetic 
terms. In the following we will show that a system with 
a Phantom and neutrinos can naturally give rise to a re- 
alization of Quintom models. Furthermore since the neu- 
trinos are the particles existing in the standard model of 
the elementary particle physics this model introduces less 
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degree of freedom in comparison to the Quintom model 
considered above with two scalar fields. 

The Lagrangian for such a coupled system is given^ by 



The equation of state for such a coupled system is 



C = + C^ + M{cj))Dv 



(14) 



where 

C^^~]^d,(pd^cp-V{cp), (15) 
with y((/)) the potential for Phantom field and 



(16) 



which is the kinetic term of the neutrino. In (|14|l M{(j})vi' 
characterizes the interaction between the Phantom and 
the neutrinos which for instance in R,ef . |2rij| is given by 

Unt = e^''^ jIlIlHH + h.c, (17) 

where /? is the coefficient characterizing the strength of 
the Phantom interaction with the neutrinos, / is the scale 
of new physics beyond the Standard Model which gener- 
ates the B — L violation, II, H are the left-handed lep- 
ton and Higgs doublets respectively. When the Higgs 
field gets a vacuum expectation value < H >^ v, the 
left-handed neutrino receives a Phantom field dependent 



f 



Majorana mass ~ e 

In general, with different V{(l)) and M{(j)) in (14) and 
(15) this two-field Quintom model of dark energy behaves 
differently. Firstly we will analyze the general feature of 
the cosmological evolution for this class of models with- 
out specifying the explicit form of V{(j)) and M{(j)). Due 
to the existence of the interaction term M{4)), the evolu- 
tion of the Phantom field is determined by the effective 
potential which is the combination of the potential V{(t>) 
andlii 



V{(j)) = nM{(j)){ 



W), 
E ' 



with n and E being the number density and energy of the 
neutrinos respectively and () indicating the thermal av- 
erage. For relativistic neutrinos, the term V^cj)) is greatly 
suppressed and the neutrinos and dark energy decouple. 
For non-relativistic neutrinos, the effective potential of 
the system is given by Ve//((/>) = V{(t)) + nM{(j)). Con- 
sequently for the equation of motion of the scalar field (/) 
it is 



(19) 



^ The studies on the int eracting Phantom dark energy with dark 
matter are given in l26l 



■U^ + V{(l))+nM{(P)' 



which can be rewritten as 



(20) 



(21) 



From the equation above one can see that during the 
radiation dominant period since Qi, is much larger than 
170, w of the coupled system will be around 1 /3 and in the 
matter dominant regime w will be close to 0. However in 
the late time when the Phantom energy dominant over 
neutrinos fi^ ^ $7^, w gradually evolves into the value 
smaller than -1. To illustrate this behavior in Fig. 4 
we plot the evolution of the equation of state w of a 

model where V — Vqc ^ "p' and M = Me ^ ""''p' . In the 
numerical calculation we have taken that A = 5.5 and 
7 = 2.5. One can see from this figure that w changes 
from above -1 to below -1 as the redshift decreases. 




Ina 



(18) FIG. 4: Plot of the w of the system as a function of the scale 
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factor Ina for V = Vbe ' "p' and M = Me ' ""'p' 



IV. CONCLUSION 

In this paper we have studied theoretically the possi- 
bility of building the two-field dark energy model with an 
equation of state across -1. In general within the frame- 
work of the general relativity and the field theory, it is 
difficult to realize the Quintom with a single scalar field 
model. Thus two fields are required. In models with 
two scalar fields there will be a lot of possibilities for the 
model building by choosing the potentials, interactions 
and the model parameters. As shown in Figs. 2 and 3 
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the physics associated with this type of models and their 
imphcations in cosmology are quite rich. In the model 
with Phantom and neutrinos, as we emphasized in the 
paper since the neutrinos are the particles existing in the 
standard model of the elementary particle physics, com- 
pared with the model with two scalar fields this model 
introduces less degree of freedom. Furthermore the neu- 
trino masses vary during the evolution of the Phantom 
field, which makes this model more interesting. Before 
concluding we should also point out that there exist pos- 
sibilities of building models of dark energy with equation 
of state across -1 beyond the field theory with minimal 



couplings to the gravity and the four dimension|2] 
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